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1 Introduction 



One of the simplest solvable one-dimensional potentials in quantum mechanics is 
the square well potential Q . It describes the confinement of a particle trapped in a 
box with infinite walls. As is well known, the energy spectrum of such a particle is 
entirely discrete and nondegenerate. 

The square well potential problem has been examined in a supersymmetric 
(SUSY) context too ||. One finds that a sequence of Hamiltonians can be gener- 
ated starting from a free-particle potential inside the well. Such a hierarchy, which 
is controlled by a family of sec 2 -like potentials, enjoys the property that its adjacent 
members are SUSY partners. 

Recently, the idea of PT symmetry in quantum mechanics has evoked a lot 
of interest ||, |], [5], [?J f|, especially in connection with its nontrivial role in non- 
Hermitian SUSY systems || [K| |ll|, [12[]. Briefly, PT-symmetric Hamiltonians are the 
ones which are invariant under joint action of parity (P: x — > —x) and time reversal 
(T: i — > — i). More importantly, such Hamiltonians are conjectured |L3|, [14], [T^, fy|, 
[17], |IJ to preserve the reality of their bound-state eigenvalues, except possibly when 
PT is spontaneously broken. SUSY methods enable one to construct non-Hermitian 
Hamiltonians with real and/or complex discrete eigenvalues by complexifying the 
underlying superpotential. Indeed this has allowed extensions of a number of exactly 
solvable models in ordinary quantum mechanics to the non-Hermitian sector || [L2|] . 

In this Letter, our primary concern is to derive a PT-analogue of the aforemen- 
tioned sec 2 -hierarchy. In this regard, the PT-symmetric version of the square well 
potential and its associated SUSY partner prove to be our natural starting point. 
Enquiries into the functioning of the PT-symmetric square well reveal [19|, |2(| that 
a critical value of the coupling parameter exists beyond which some energy levels 
appear in complex-conjugate pairs. This signals an onset of a PT-spontaneously 
broken phase. However, below the critical value, PT symmetry is unbroken with 
all the energy levels remaining real. As we plan to show, the existence of a PT- 
symmetry broken phase is responsible for some new features in the SUSY extension 
to non-Hermitian Hamiltonians. 
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2 PT- Symmetric Square Well 

A PT-symmetric square well potential on the interval (—1, 1) can be defined in the 



manner 
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V R {x) = -iZ, V L (x) = iZ, (1) 

where Z is the coupling strength and L (R) denotes the region — l<x<0(0<x< 
1). The wave functions, at the end points of the interval, are enforced to be vanishing 
as is normally the case with the real square well problem with impenetrable walls. 

For real eigenvalues E n (n = 0, 1, 2, . . . ), the Schrodinger equation for ([[]) is 
equivalent to the pair 

<R = «WV*, €L = <VnL, (2) 

where 



n 2 n = -E n - iZ = (s n - it n ) 2 , s n , t n e E. (3) 



Solving @ gives 



t n = ^= (E n + sjEl + Z^j , s n = -| (E n + sjEl + Z^j ~ 1/2 . (4) 

The solutions of (0) fulfilling the conditions t/VirXI) = i>nL{— 1) = may be 
written as 

?pnR = C nR sinh[/t n (l - x)}, ip nL = (7 nL sinh[/«*(l + x)]. (5) 

To establish a link between the (complex) constants C„r and C„l, we see that the 
continuity of the wave function and its derivative at x = imposes the conditions 

K n coth K n + K* n coth K* n = 0, (6) 
C nR sinh < 



C nL sinh k„ 

Now, if we require ip n n and t/Vl to be also PT-symmetric near the origin, that is 



(7) 



^nR(O) = tpnUO) = a n , d x ip nR (0) = d x ?p nL (0) = ip n , (8) 
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where the parameters a n , (3 n G R, then use of (|[) enables us to rewrite the eigen- 
functions (||) as 

= . sinh[/t n (l - x)}, ip nL = , " n sinh[K* (1 + ar)]. (9) 
smh/c n sinhK* 

It is to be noted that the condition @, which is independent of C„r and C n L, 

can also be put in the form 

s n sinh2s n + t n sin2t n = 0, (10) 

where we have used (|3|). Setting s n sinh 2s n = 2 sinh 2 S n and t n = |T n , we can recast 
©as 

4 sinh 2 S n = —TiT n sin nT n . (11) 

Equation (|Tl~l) shows that the points (X„, 5 n ) of the T-S plane corresponding to the 
eigenvalues E n belong to the Z-independent curve 

S = X(T) = arcsinh \\f— ttT sin 7rT, (12) 

where 2m — 1 < T < 2m, m = 1, 2, . . . . On the other hand, Eq. (^) may be exploited 
to obtain sinh 2 S n = -^f - sinh which in turn exposes the Z-dependence of S: 



/ Z 2Z 

S = Y(Z,T) = arcsinh W — sinh— . (13) 

The pair of equations (|T2|) and (fl3f) imply that the points (T n , S n ), and in con- 
sequence the energies E n , are at the intersections of the two curves 5* = X(T) and 
S = Y(Z, T). As shown in jg, |2T| , if Z is below the critical threshold Z { Q CT ' lt) ~ 4.48, 



there are two real eigenvalues in every interval 2v — 1 < T < 2u, where u = 1, 
2, . . . . At Z = ZQ Crit ' ) , the two lowest real eigenvalues Eq and i?i merge and for 
Z > Zq'"^ move into the complex plane, where they become complex conjugate. 
The corresponding eigenfunctions ipo an d ipi then cease to be PT-symmetric. More 
generally, there exists a naturally-ordered increasing sequence of critical couplings 

^(cnt) < . . . a ^ which a pair of real eigenvalues (i?2i/> E 2u +i) 
merges and moves into the complex plane. 

In Sees. 2-5, we shall consider the unbroken PT-symmetry case 
leaving the discussion of the broken PT-symmetry one (Z > 4 Crit) ) for Sec. 6. 



3 SUSY Partner in the Unbroken PT-Symmetry 

Case 



With the above preliminaries on the PT-symmetric square well potential, we proceed 
to the construction of its SUSY partner in the Z < Zq 1 ' case. Adopting the 



notations of [12], we consider the following factorization scheme 



# (±) = ~ + V^\x) -Eq= (AA,AA), (14) 
ax z 

corresponding to an arbitrary factorization energy E = Eq. In ([14]) , are the 
SUSY partner potentials and the operators A, A may be defined in terms of a 
superpotential W(x) as 

A = 4- + W(x), A = -4- + W(x). (15) 
ax ax 

Inserting (15) in ([14]) gives 

V {±) = W 2 tW + E . (16) 

Identifying y( + ) with the square well potential (|l|), that is, taking = —iZ 
and = iZ, we obtain for W(x) the following differential equations 

Wl - = 4, Wl -Wl = nf. (17) 



Solving (p?7f ) we get 



Wr(x) = — K tanh[K (x — xr)], Wi,(x) = — tanhf^^x + xl)], (18) 

where xr and xl are two integration constants. 

From (|18D it results that the partner potential to ([D is given by 

V^~\x) = -iZ + 2kI cosech 2 [/t (l - x)}, U L ( ~ } (x) = iZ + 2k* 2 cosech 2 [K*(l + x)}. 

(19) 

The partner Hamiltonian is therefore PT-symmetric. Note that in (0), xr 
and xl have been chosen as xr = 1 — m/ (2« ) and xl = 1 — ht/ (2/cJ) to ensure that 



\x) blows up at the end points x — — 1 and x = +1, which is as it should be. 
Observe that the superpotential W (x) also blows up at these points: 



W R (x) = k coth[re (l - x)] f W L (x) = -k* coth[«JJ(l + x)}. (20) 

To obtain the eigenfunctions associated with H^~' we first remark that the 
ground state eigenfunction ipQ + \x) of H^ + \ given by (§) when n = 0, is annihi- 
lated by the operator A: 



V'or 0*0 = S ^_ K \~^ + K ocoth[/t (l - x)) 1 sinh[K (l - x)] = 0, 



(21) 

where the superscript (+) is appended to the eigenfunction and the coefficient pa- 
rameter to signify that we are dealing with the component. A similar result 
like (H) holds for [■£ + W L (x)] ^\x). 

Exploiting then the intertwining character of SUSY, the eigenfunctions ip^\x), 
n — 0, 1, 2, . . . , of are obtained by application of A on subject to the 
preservation of the boundary and continuity conditions: 

^r } (1) = 0, ^-)(-l) = 0, (22) 
Vfe } (0) = V& } (0), d x ^(0)=d x ^\0). (23) 

We get in this way 

V>£r 0*0 = C { ~2 g^"^ 1 sinh[/c n+ i(l - x)} 

x {-K n+l coth[/c n+ i(l - x)\ + kq coth[« (l - x)}}, 

1>ti(x) = ^-|^sinh[< +1 (l + x)] 

x {< +1 coth[< +1 (l + x)] - K* Q COth[K*(l + x)]}, (24) 

where and C^J are (complex) constants. These eigenfunctions satisfy 

the boundary conditions (p2|) because of the limiting relations of the type 
lim a; _ > i{sinh[/c n+ i(l - x)\/ sinh[K (l - x)}} = k u+1 /k . 
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On the other hand, because of the continuity conditions (|23| ) it turns out that 

CnR = ^nh = Cn \ (25) 

along with 

K n+1 ~ K n+1 = K ~ K Q = — 2iZ, (26) 

where we have used ([|). 

As previously, if we require ip^ to be also PT-symmetric near the origin and so 
denote by and the respective values of i}j(~~'(x) and d x ip!^\x) at x = 0, 
where a^~\(3^ E M, we obtain, on account of (|25j), the following forms of the 
eigenfunctions of H^~': 



a4 ) sinh[K n+ i(l - x)} 



sinh K n+ i(K n+ i coth « n+ i — k coth k ) 

x {n n+l coth[/c n+1 (l - x)] - k coth[/c (l - x)]}, 



{ _) _ o4 )sinh[< +1 (l + x)] 



sinh k* +1 (k* +1 coth — Kq coth Kg) 
x coth[< +1 (l + x)] - k* cothK(l + x)]}, (27) 

where c^"- 1 = C^~- ) a^ f \(— /t n+ i cothre^+i + n coth.K, ). 

We remark that the above eigenfunctions of H^~\ defined by ([14]) and (0), have 
SUSY related eigenvalues 

= E n+l = E n+1 - Eq = K 2 Q - K* +1 . (28) 



Notice that in (p8|), the coupling strength Z only appears implicitly in the «'s. 
Equation ( p8|) reflects a typical unbroken SUSY feature: pairing of the eigenvalues 
of the partner Hamiltonians with the ground state nondegenerate for n = 0, as 
shown by (plj) for ip^(x) and a similar equation for ip^ (x). 

4 Z -> Limit 

At this stage, it is instructive to look into the Z — > limit. In this limit, Eq. (H) 
becomes k 2 = —E n = —t\. From the continuity condition (|H]), which is now 



t n sin2t n = 0, we are led to 



(n + l)|, K = (n + 1) 2 ^. (29) 
7 



For the odd values of n for which sint n = 0, the other continuity condition ([?]) 
becomes useless because its right-hand side is indeterminate. Going back to the 
square well eigenf unctions (H), which are now 



^ = -i^sin[t n (l-x)], 



»(+) 



(+) 

riL 



iCit } sin[t n (l + x) 



(30) 



it is however straightforward to see that the continuity conditions yield two solutions 



C 



c 



(+) 

nR 

(+) 
nR 



_ r (+) _ W+) 

M+) _ r (+) 



cost n = 0, 
sint n = 0, 



(31) 



according to whether n is even or odd. As a consequence, the eigenfunctions can be 
written as 



4% 



-id +) < 



-'2;/ 



-ireos 



7T 



id +) 



2v+V 



(2u + l)-x 
sm[(v + 1)ttx], 



(32) 



where we do not have to distinguish between the intervals (—1, 0) and (0, 1) anymore. 
The forms (|32|) are in conformity with the known results for the real square well Jl| . 

Let us now consider the SUSY partner as Z — > 0. Since kq = — ii = — i|, 
the superpotentials in (p0|) become W-^^x) = | tan (fa;). As a result, the partner 
potentials in (|T§]), too, acquire the common form = ^ sec 2 (fa;) that coincides 
with the SUSY partner of the real square well first obtained in ||. 

For the partner eigenfunctions we obtain, in the Z — > limit, 



C^Cfik £{-\y +l {(2u + 2) «*[(,/ + l)nx] 
+ tan f^ 2 ^ sin[(z/ + l)7rx]|, 
cU 1 cW 2 i^(-l)^ 1 {(2z. + 3)sin 



(2u + 3)-x 



— tan ( —x ) cos 

.2 



7T 



(2z/ + 3)-x 



From the results of ||, we know that the real square well V\(x) - 
generates a whole family of sec 2 -like potentials with increasing strengths, 



7T 2 / 7T 

V m {x) = V%(x) + —m{m — 1) sec 2 I — x 



m = 1, 2, 3, ... , 



(33) 
(34) 



corresponding to a hierarchy of Hamiltonians, whose adjacent members are SUSY 
partners. Having found the PT-symmetric analogue of the second member of the 
family, V^ix) = (x), we may now try to build counterparts of the other members, 
V3, V4, .... Such a construction is outlined in the next section. 

5 SUSY Hierarchy in the Unbroken PT- 
Symmetry Case 

Let us define a hierarchy of partner Hamiltonians H m , m = 1, 2, . . . , whose first 
member Hi coincides with that of the PT-symmetric square well. According to this 
description, we have the following set of SUSY partners 



— ~~i^2 ^ / m^^i x ) ~ E m Q — H m — E m fl — A m A m , 
d 2 

Htf = -— + Vi-\x)-E mfi = H m+1 -E mfl = A m A m , (35) 
where V^\x) = V m (x), V^\x) = V m+ x(x), m = 1, 2, . . . . For m = 1, Vir and 



Vil are given by ([!]), while for m = 2, U 2 r and V 2 l are given by fljjp . As usual, the 
operators A m and A m can be written in terms of the superpotentials W m and are 



A ™ = £ + W m (x), A m = -£ + W m (x), m = 1, 2, . . . . For m = 1, W m and W lh 



We denote the eigenvalues and eigenfunctions of H m by E m ^ n and Vv«> respec- 
tively. These satisfy the SUSY properties 



are given by (POD- In terms of W m , the partner potentials Vj^ 1 ' read 



V {±) = W 2 =fW' +E m n . 



(36) 





m-l,n+ljR„L, 



(37) 



along with the conditions 



(VwOr(°) 



9 x (^m,n)R(0) 



5x(-0m,n)L(O) = i/3 m , n , 



(38) 
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where a mj „, (3 m ^ n G E. Up to now we have determined the energy eigenvalues for 
m = 1 and 2, 

E hn = E n = -iZ - K 2 n , E 2: n = -iZ - K 2 n+1 , (39) 
with associated eigenfunctions 

Oi,n)it = (Ci ifl ) R sinh[K n (l - x)], 

C0i,n)L = (Ci )fl ) L sinh[/c*(l + 

(V ; 2,n)R = (C , 2 in ) R (C , i in+ i) R sinh[/t ri+1 (l - a;)] 

x {—K n +i coth[/t n+1 (l — x)] + kq coth[/«o(l — ^)]}, 

(■02,n)L = (C2,n)L(Ci, n+ i) L sinh[K* +1 (l + x )] 

x {< +1 coth[< +1 (l + a;)]-K*cothK(l + a;)]}. (40) 

In ( [40|) we have defined 



sinh K n ' sinh /c* 

(C 2 ,n)R = (C 2 ,n)L = C2,n = 7 n~ i 7u V' 

ai, n +i(-Kn+i coth K n+ i + K coth «o) 
To construct the third member of the hierarchy we observe that Wz(x) can be 
easily determined from the property A 2 ip2,o = that yields W 2 (x) = — ip' 2 o/V^o- 
Explicit calculations give for W 2 the forms 



W 2>R (x) = -«ocoth[«o(l-a:)] + r . r r 

«i coth[«i(l — x)J — /to coth[/«o(l — £)J 

^ 2 ,l(x) = -W* R (-x). (42) 

For the potential V 3 (x), our results turn out to be 

2 2 k\ cosech 2 [/ti(l — x)\ — Kq cosech 2 [K (l — x)] 
V 3>R {x) = -iZ - 2{k x - k ) { ^ CQth[Ki(1 _ x)] _ Kq coth[Ko(1 _ ^jp ' 

^3,l(x) = V£ R (-x), (43) 
while for the eigenfunctions we get 

C03,ti)r = (C3 in ) R (C2,n+i)R(Ci, n+ 2)Rsinh[K ri+ 2(l - z)] j/4+2 - 
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n L 



K\ coth[/ti(l — x)] — K coth[/t (l — x)] 

(C , 3,n)L(C , 2,n+l)L(Ci in+2 )L Smh[< +2 (1 + x)]\ < 2 +2 - K* Q 2 



*2 



<2 < +2 coth[< +2 (l + x)] - coth[/cS(l + x)\ 



(44) 



/?* coth[/«i(l + x)\ — Kq coth[/«o(l + ^)] 

where E 3j „ = -iZ - «^ +2 and (^ 3 ,n)R(l) = = 0. 

It remains now to impose the continuity conditions (|38|) on the eigenfunctions 
HI). While the first one leads to 



{Cz,n)-R — {C-i,n)h 
"3,n 



a 



3,71 



v n+2 



(45) 



«2,n+l I K l COth «i — K COth Kq K n +2 COth K n+2 — ^0 COth Kq 

the second one amounts to an identity when (|4"I) and (33) are taken into account. 

We have thus obtained explicit forms of the first three members of the SUSY 
hierarchy for the PT-symmetric square well potential. It is clear that by applying 
similar techniques, formulas for other members may be similarly constructed. 

It is easy to check that in the Z — > limit, all the results obtained in this section 
go over to those for the third member of the real square well hierarchy. Equations 
( f42"D and (f43|), for instance, yield W 2 ^{ x ) = W 2 ,-l{x) = W 2 (x) = vrtan^fx) and 
^,r(^) = ^3,l(^) = V 3 (x) = |7r 2 sec 2 (fx), in conformity with Eq. (|34]) . 

It is useful to stress here that, in the same limit, the eigenfunctions of the first 
three members in the hierarchy, namely (Q), (|27|), and (|44|), turn out to be propor- 
tional to Gegenbauer polynomials: 



: 



-\Cxjn COS ( ^X ) C£ ] 



71 

Sill I — X 



-i7rC 2 ,„Ci >n+ i cos 2 ( ^x ) Cf> 



7T 

sin I — x 



-2i7r 2 C 3 ,nC 2in+1 C lin+2 cos 3 f |x) Cf 



7T 

sin I — x 

2 



.(46) 



In (|32|) and ( |33l) we had already furnished the limiting forms [Z — > 0) of the first 
two members in the hierarchy. Note that to get to the representations ([46]) we used 
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the definition of the Gegenbauer polynomial C^ 1 ** (cos (/>) = [sin(n + 1)0]/ sin0 and 
considered the general recurrence relation 



d 7T / 7T 

- — h (m — 1)— tan — x 

dx v ; 2 V2 



7T 

sin 1 — x 
2 



n(m - 1) cos m (|ac) 



71 

SU) | — X 



(47) 
(48) 



which can be easily obtained from known properties of Gegenbauer polynomi- 
als For nonvanishing values of Z, it can be shown that the eigenfunctions 
of the PT-symmetric square well and of the next two members in the hierarchy, 
namely ( [TOD and (|44D , can be rewritten in terms of Gegenbauer functions of the 
type C ( ( ^ WKo) _ p {cosh[/€ (l - x)}} or C ( ( ^ /k5) _ p {cos1i[k*(1 + x)]}, where p,qEE. 

6 SUSY Hierarchies in the Broken PT-Symmetry 

Case 



Let us now consider coupling strengths Z for which PT symmetry is spontaneously 
broken and assume first that Z lies between the first two critical values, Z^ 1 ^ rs 



4.48 < Z < Z\ 



(crit) 



12.80 



The PT-symmetric square well has then a single 



pair of complex-conjugate eigenvalues E = eo — ieo, E\ = Eq = eo + ieo (where 
e > 0), and an infinite sequence of real eigenvalues E n , n = 2, 3, . . . . 
The corresponding Schrodinger equation is equivalent to 



(49) 



where 



Po 



A 



o- 



Pi = A 0) p 

(Tl = AC , 



n = 2,3, 



n = 2, 3, 



(50) 



and 



-e + ie - iZ, A = -e - ie - iZ, 



-E n — iZ, n = 2, 3, . . . . 



(51) 
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The solutions of P9|), vanishing at the end points of the interval (— 1,+1) and 
satisfying the continuity conditions of ip n and ip' n at x = 0, are given by 



ipnR = C nR sinh[p n (l -x)], 



C nL sinh[cr n (l + x)], 



(52) 



where 



p n coth p n + a n coth a n = 0, 




sinh a. 



sinh p. 



(53) 



Choosing them real at x = 0, we obtain 



a. 



•n 



a 



(54) 



sinh p. 



sinh a. 



n 



With this choice, the eigenf unctions ip n , n — 2, 3, . . . , corresponding to the real 
eigenvalues E n , n = 2, 3, . . . , are PT-symmetric as in the unbroken PT-symmetry 
case, but this does not hold true for ip , ipi, corresponding to the complex-conjugate 
eigenvalues i^o, E\. 

To construct a SUSY hierarchy of partner Hamiltonians H m , m = 1, 2, . . . , 
whose first member coincides with that of the PT-symmetric square well, we have 
now various possibilities at our disposal. At each stage, we may indeed eliminate 
either the lowest real eigenvalue or one of the two complex eigenvalues as long as 
there remains some. As we shall proceed to show, this leads to a rich diversity 
of hierarchies: a PT-symmetric hierarchy with spontaneous symmetry breaking, 
two non-PT-symmetric hierarchies with partial PT-symmetry restoration, as well as 
various mixed-type hierarchies. 

Let us start with the first one to be distinguished by an upper index (1). By 
successively eliminating the lowest-lying real eigenvalues E^l = E2, E^l = E3, . . . , 
while keeping the two complex ones E^q = E , E^l = E\ at every stage, we obtain 
PT-symmetric potentials V^'(x), m = 1, 2, 3, . . . , with a pair of complex-conjugate 
eigenvalues E^ = E , E„\\ = E\, and an infinite sequence of real eigenvalues 
E^' 2 — E m+ i, Ej^\ = E m+2 , .... The corresponding eigenf unctions ip$ n (x), n = 0, 
1, 2, . . . , turn out to be PT-symmetric only starting from n = 2. 
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This is illustrated with the second and third members of the hierarchy, which can 
be obtained from the results of Sec. 5 by making some appropriate substitutions: 
K o Pi-, K i P3? K o — * °"2! K i - * °3 for Wi~\ V^, W^~\ and V^, and the same 
together with K\ — * p , k 2 — > pi, — > o"o ? K 2 — * °i or K 2 - Po, K 3 — ► Pi) K 2 — * a 0) 
«3 — > o"i for an d V4*n) respectively. The corresponding potentials, for instance, 
are given by 

V^r(x) = — iZ + 2k 2 , cosech 2 [/t 2 (l — #)], 

= ^'(-s), (55) 



and 



(1) 2 2 /t|cosech 2 [/t 3 (l -x)] -^cosech 2 [/t 2 (l-a;)] 



{ks coth[/t 3 (l — x)] — k, 2 coth[/t 2 (l — x)}} 2 ' 
= (56) 

We conclude that all the numbers of this first hierarchy have properties very similar 
to those of the PT-symmetric square well for the chosen Z value, namely all of them 
are PT-symmetric but exhibit spontaneously broken PT symmetry with a single 
pair of complex-conjugate eigenvalues. 

Let us now consider the second type of hierarchies, which are obtained by suc- 
cessively eliminating the two complex eigenvalues, then the real ones in increasing 
energy order. Since we may choose to eliminate first Eq, then Ei, or the reverse, 
there are two different hierarchies to be referred to by an upper index (2) or (3), 
respectively. 

To obtain results for the former from those of Sec. 5, it is enough to perform the 
substitutions K n — > p n , «* — > a n . So we get for instance 

^2,r( x ) = — iZ + 2/€q cosech 2 [/t (l — a;)], 

V?i(x) = iZ + 2A* 2 cosech 2 [A*(l + x)], (57) 



and 



(2) _ _- 7 _ 9 / X 2_ 2 ^ A 2 cosech 2 [A (l - x)] - /t 2 cosech 2 [/t (l - a;)] 
V 3jR {X) - iZ Z{A K ) {XoCoth[Xo{1 _ x)] _ KoCoth[Ko{1 _ x)]} 2 > 

V$(x) = V™*(-x), (58) 
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while the remaining potentials Vjf\x), n = 4, 5, are PT-symmetric. The 

(2) (2) 

potential V 2 (x) and the eigenfunctions of the corresponding Hamiltonians H 2 
are not PT-symmetric, which is not surprising ||, [14|, [l^, |18|, [17], |T6[| since 



2 



(2) (2) 
has a single complex eigenvalue E 2 q = E\ in addition to real ones, E 2 ^ = E n+ \, 

(2) 

n = 1, 2, . . . . Strangely enough, the next potential V 3 is PT-symmetric and the 

(2) 

corresponding spectrum E^ = E n+2 , n — 0, 1, 2, . . . , is entirely real, but the 

(2) 

eigenfunctions ^^(x), with no definite symmetry under PT, do not even differ by 
a simple phase factor from ip^(-x), as it is normally the case for eigenfunctions 
of PT-symmetric Hamiltonians [|HJ . This distinctive feature may be traced back to 
the finite discontinuity of the PT-symmetric square well and of its SUSY partners 
at x = 0, which has led us to ensure smoothness of the eigenfunctions by imposing 
continuity conditions externally by hand. As explicitly shown here, such a procedure 
does not guarantee PT symmetry of the eigenfunctions at all levels. In contrast, the 
results of JT5J were based on the tacit assumption of smoothness of the potential over 



the entire real line and may therefore not be compared with the present situation. 

The other second-type hierarchy only differs from the first one by the interchange 
of E Q and E\ or k , an d -^o, ^o- As a result, we get 

V 2 {3 \x) = V 2 (2) *(-x), Vf } (x) = Vj 2 \x). (59) 

Hence all the corresponding potentials of the two hierarchies coincide, but for the 
second ones, which are related through PT symmetry. It can also be shown that 
from m = 2 onwards, the eigenfunctions of corresponding members Hlff of the 
two hierarchies are also related through PT symmetry: ^$ n (x) = ijj^*(—x), m = 2, 

3, . . . , n = 0, 1, 2, We therefore conclude that these second-type hierarchies 

are non-PT-symmetric ones, but exhibit partial PT-symmetry restoration. 

It should now be clear that apart from the hierarchies considered so far, there 
also exist a lot of mixed-type ones, which differ from one another and from the 
previous ones in the step (resp. steps) where one of the two (resp. both) complex 
eigenvalues is (resp. are) eliminated. 
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Finally, if the coupling strength Z lies in another interval \Z^\ , Z^^J, such 
that v e 2, 3, . . ., and there therefore exist v pairs of complex-conjugate eigenvalues 
(E , E 1 = £*), (E 2 , E 3 = £*), (E 2v - 2 , E^-i = E* v _ 2 ) in addition to the 
real ones, SUSY hierarchies of partner Hamiltonians can be constructed in many 
different ways. In particular, by eliminating the complex eigenvalues, it is possible 
to partially restore PT symmetry after 2is steps. 

7 Conclusion 

In the present Letter, we have established the SUSY connection of the PT-symmetric 
square well both in the unbroken and broken PT-symmetry cases. In the former 
we have availed ourselves of this to derive a PT-symmetric analogue of the sec 2 - 
hierarchy. 

In this respect, the PT-symmetric world proves more intricate than the Her- 
mitian one: not only has one to resort to numerical calculations to determine the 
eigenvalues, but also the potentials and the eigenfunctions get more and more com- 
plicated when going to successive members of the hierarchy in contrast to what 
happens for the real square well. 

Another intricacy of the PT-symmetric square well problem, namely the exis- 
tence of complex eigenvalues for a coupling strength above the critical threshold 
ZQ Tlt \ leads to a new feature in the SUSY extension to non-Hermitian Hamilto- 
nians: the existence of a rich diversity of SUSY hierarchies, including some with 
PT-symmetry breaking and some with partial PT-symmetry restoration. 
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